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.3  abstract  -phe  raCi0  Gf  o / T^  for  mixed  and  free  convection  has  been  examined  as  a 
function  of  heat  flux,  or  t h/b  % of  inversion  height,  and  of  the  effects  of  water  vapor 
on  buoyancy.  For  the  stability  range  -Z/L>0.1,  the  free  convection  similarity  relation 
o /Tv=-0.95(-Z/L)-l/3  is  accurate  up  to  at  least  60%  of  the  inversion  height.  For  t lie 


lower  heights  where  forced  convection  is  important  the  relation:  o /T.=-0.95(0.0>49- 

i / l t * 


Z/h)“^/3  provides  a stability  dependent  correction  for  the  heat  flux.  The  effect  of 
water  vapor  on  the  buoyant  processes  can  be  incorporated  by  substituting  virtual  temp- 
erature for  temperature  in  the  relations  above.  These  relations  also  determine 
implicitly  the  buoyancy  flux  at  the  surface.  This  is  a good  approximation  of  the  sen- 
sible iieat  i 1 ux  as  long  as  the  Bowen  ratio  is  larger  than  one. 

The  horizontal  uniformity  of  minimum  temperature  during  quiescent  intervals  for  unstabl 
convection  is  shown  to  be  on  the  order  of  a few  hundreds  of  a degree  celsuis  over  a 
fifty  meter  array  for  periods  on  the  order  of  fifteen  minutes.  One  interpretation  of 
this  characteristic  is  that  the  minimum  temperature  during  the  quiescent  interval  is 
characteristic  of  the  potential  temperature  at  planetary  boundary  layer  heights  .it 
which  the  potential  temperature  is  constant  with  height.  Examination  of  the  quiescent 
temperature  during  unstable  and  stable  periods  of  the  same  day,  shows  that  the 
spatially  uniform  quiescent  temperatures  go  from  cold  during  unstable  conditions  to  war 
during  stable  conditions.  The  latter  is  consistent  witli  the  hypothesis  that  elements 
from  height  ranges  of  uniform  potential  temperature  are  adverted  into  the  surface  layer 
by  mechanical  mixing  during  these  stable  conditions. 
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Spatial  Temperature 


FIELD  PROGRAMS 


Three  sets  of  data  were  either  partially  or  completely  acquired  and  analyzed 
under  this  grant-  The  major  field  program  funded  was  the  Minnesoto 
experiment,  conducted  in  the  fall  of  1973  in  northwestern  Minnesoto.  (Izumi 
and  Caughey.  1976)  The  University  of  Washington  group  joined  the  experiment 
as  a guest  of  the  United  States  Air  Force  Cambridge  Research  Laboratories' 
Boundary  Layer  Research  Croup,  under  the  leadership  of  Dr.  Duane  Haugen. 
This  cooperative  mode  of  experiment  was  chosen  due  to  the  small  budget 
required  for  the  field  progam  and  the  anticipation  that  a major  portion  of 
the  analysis  required  to  provide  the  parameters  that  characterize  the 
planetary  surface  and  boundary  layer  would  be  a by  product  of  the  AFCRL 
experiment.  Another  reason  for  utilizing  the  AFCRL  facilities  is  that  the 
site  was  chosen  for  horizontal  uniformity.  In  many  respects,  the 
anticipated  advantages  were  not  realized.  Our  participation  in  the 
experiment  was  delayed  until  the  other  groups  had  obtained  their  data.  The 
AFCRL  group  lost  two  balloons  due  to  meteorological  phenomena,  causing  them 
to  terminate  the  field  program  earlier  than  expected.  The  anticipated  data 
analysis  benefits  were  not  realized  as  the  efforts  were  devoted  to  the  major 
joint  experiment  with  the  British  Meteorological  Office.  Since  problems 
were  anticipated,  we  obtained  real  time  output  of  15  minute  average  profile 
and  flux  quantities.  Due  to  the  scatter  of  the  real  time  field  estimates  of 
the  boundary  layer  parameters  from  the  measured  fluxes,  we  obtained 
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estimates  from  the  profiles  of  temperature  and  wind  by  analyzing  the  profile 
wind  and  temperature  data  using  the  flux-gradient  relations.  These 
quantities  became  the  basis  for  our  estimates  of  the  boundary  layer 
parameters  after  eliminating  data  from  levels  on  the  tower  with  faulty 
sensors.  An  added  benefit  of  the  program  was  obtaining  the  data  from  the 
balloon  which  provided  temperature  statistics  throught  the  planetary 
boundary  layer-  Measuring  the  height  of  the  inversion  and  it's  effect  on 
the  parameters  was  an  important  factor  in  the  experiment  since  this  is  our 
only  data  set  over  land  where  the  inversion  height  is  available.  However, 
the  data  from  this  experiment  taken  at  the  heigher  levels  is  anomolous  which 
led  to  a more  thorough  investigation  of  the  y /T*  relation. 

The  site  characteristics  were  not  as  ideal  as  expected  since  the  extreme 
flatness  of  the  area  made  slight  indentations  form  muddy  areas  which 
remained  an  anomolous  source  of  water  vapor  flux  surrounded  by  dry  soil. 

This  may  be  responsible  for  some  of  the  puzzling  results  from  one  or  two 
sensors.  However,  the  uniformity  in  roughness,  with  proper  fetch,  provided 
a good  setting  to  verify  the  uniformity  of  potential  temperature  during 
quiescent  intervals  between  the  Douyant  convective  elements. 

J 

The  puzzling  aT/T*  versus  Z/L  relation  for  the  Minnesoto  balloon  data 
warranted  furthur  investigation  especially  since  this  relation  is  one  of  the 
most  stable  and  site  independent  encountered  prior  to  this  experiment. 
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Since  the  data  in  question  was  obtained  at  heights  that  are  an  appreciable 
fraction  of  the  inversion  height,  as  opposed  to  that  of  the  tower  based 
experiments  data  covering  the  range  of  planetary  boundary  layer  heights  was 
obtained  and  analyzed.  ( Hiester,  T.  R..  1977  ) 

An  objective  of  the  program  was  to  determine  to  what  degree  of  accuracy  the 
quiescent  temperatures  observed  during  highly  unstable  temperatures 
represented  the  potential  temperature  observed  at  a significant  fraction  of 
the  inversion  height.  The  Buffalo  aircraft  of  NCAR  was  scheduled  to  fly 
during  our  phase  of  the  experiment  to  permit  direct  comparison  of  the 
surface  layer  quiescent  temperature  with  that  measured  over  time  and  space 
from  the  aircraft-  Unfortunately,  structual  problems  with  the  aircraft 
cancelled  this  part  of  the  program.  However,  the  space  and  time  variability 
within  the  array  has  been  used  to  support  the  hypothesis  that  the  parcels 
with  quiescent  temperatures  originate  significantly  above  the  surface  layer. 

A primary  objective  of  the  Minnesota  temperature  array  analysis  was  to 
determine  the  translation  velocity  of  convective  plumes.  Ideally,  the 
direction  of  the  wind  at  the  time  of  passage  of  a plume  should  be  used  for 
the  translation  velocity.  However,  because  the  time  series  of  wind  speed 
and  direction  were  not  available,  it  was  necessary  to  make  some  assumptions 
concerning  the  translation  direction  of  the  plumes.  As  reported  previously 
( our  and/or  other  notes  and  papers  ).  assuming  that  the  plumes  translate  in 
the  direction  of  the  mean  wind  leads  to  substantially  greater  translation 
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velocities  than  if  it  is  assumed  that  the  plumes  translate  in  a direction 
perpendicular  to  their  orientation.  ( Their  orientation  is  determined  by 
fitting  a plane  wave  to  the  sharp  temperature  discontinuity  at  the  upwind 
edge  using  three  sensors  in  a triangular  array.  ) Since  there  is  no  clear 
reason  to  select  one  method  of  determining  translation  direction  over  the 
other,  the  translation  velocity  results  remain  inconclusive. 

The  second  objective  of  the  analysis  was  to  determine  the  planform  pattern 
of  the  plumes  in  a horizontal  plane.  However,  to  do  so  again  requires  the 
knowlege  of  the  translation  direction,  in  addition  to  the  translation 
velocity.  It  was  therefore  impossible  to  determine  the  planform  pattern  of 
the  convective  plumes.  The  inconclusive  results  led  to  a second  experiment 
which  took  place  at  the  Roulder  Atmospheric  Observatory  in  April  1978- 
During  this  experiment,  both  wind  speed  and  direction  were  measured  along 
with  the  array  temperature  using  the  same  temperature  system  as  constructed 
for  the  Minnesota  experiment.  The  measurements  were  made  by  J.  Tillman  and 
J-  Wilczak.  Preliminary  results  indicate  that  conclusive  determination  of 
both  the  translation  velocity  and  planform  structure  of  convective  plumes 
will  be  possible  from  the  data  set  and  the  results  will  be  reported  in  a 
Master  of  Science  thesis  being  prepared  by  James  Wilczak. 
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The  layout  of  the  temperature  sensor  array  constructed  for  this  program  is 
given  in  Figure  1.  The  majority  of  the  sensors  were  located  at  a height  of 
four  meters  along  a south  to  north  and  a west  to  east  line  from  the  center 
of  the  array.  The  center  was  located  50  meters  east  of  the  AFCRI,  tower  and 
our  sensors  were  located  at  heights  of  2,4.8.16  and  32  meters  on  the  AFCRL 
tower.  At  the  center  of  the  array,  sensors  were  located  at  heights  of  0-5. 
I-  2.  and  4 meters.  The  vertical  supporting  poles  were  7.5  cm  in  diamater 
and  the  sensors  were  located  one  meter  from  the  poles  to  minimize  the 
interception  of  hot  plumes  from  the  poles.  The  sensors  are  25  micron 
platinum  wire  in  a three  wire  AC  synchronous  modulation-demodulation 
resistance  thermometer  system.  The  basic  system  exhibits  stabilities  of 
better  than  0.01  degrees  over  long  periods  exclusive  of  any  stretching  of 
the  platinum  wire.  A more  complete  description  is  given  by  Katsaros,  et  al. 
1977. 

Calibration  of  the  complete  system,  exclusive  of  the  sensors  was 
accomplished  in  the  field  by  substituting  five  precision  fixed  resistors  in 
each  channel  after  the  completion  of  the  field  experiment.  The  calibration 
data  was  written  on  tape  in  exactly  the  same  manner  as  the  data  were 
collected.  Sensors  were  calibrated  to  within  a range  of  3 degrees  absolute 
using  a hath  of  mineral  oil  and  a calibration  adjustment  was  made  using  the 
quartz  crystal  thermometers  of  AFCRL  as  a reference  over  a 15  minute  period 
during  the  field  program  to  improve  the  accuracy.  The  adjustment  consisted 
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of  assuming  that  the  electronic  system  was  accurate  ot  0.01  degrees  and  that 
any  errors  were  due  to  stretching  of  the  platinum  wire.  With  this 
assumption.  a single  parameter  adjustment  was  made  to  bring  the  mean 
temperature  of  all  array  sensors  into  agreement  with  the  quartz  thermometer 
at  the  same  height.  The  data  of  Table  1 show  the  reduction  after  this 
adjustment  where  the  AFCRL  and  UW  15  minute  mean  temperatures  are  equivalent 
by  definition.  Note  that  the  standard  deviations  of  the  sensors  in  the 
array  agree  to  within  0.045  degrees  even  though  the  calibration  was  a single 
parameter  change  rather  than  a bias  and  a gain  adjustment-  Subsequent 
checks  of  the  sensors,  after  the  elimination  of  several  noisy  sensors, 
produced  maximum  differences  of  0.12  degrees  between  any  of  the  array 
sensors  and  the  quartz  sensors  on  the  tower  except  below  2 meters  where  the 
AFCRL  tower  base  and  it's  surrounding  area  are  unrepresentative  of  the 
surface-  Generally,  the  differences  were  on  the  order  of  a few  hundredths 
of  a degree.  The  subsequent  run.  shown  in  Table  2.  2 hours  and  12  minutes 
later,  indicates  a maximum  difference  of  0-05  degrees. 

Filtering  of  the  signal  was  accomplished  by  three  pole  filters  whose 
response  was  down  by  50%  at  10  Hz.  The  ptiase  shift  was  controlled  so  that 
no  signal  channel  differed  in  phase  shift  from  any  other  more  than  10 
degrees  at  10  Hz.  Data  collection  was  accomplished  by  a Raytheon  704 
computer  writing  standard  7 track  tapes.  Data  sampling  rates  were  program 
selectable  and  there  was  no  limit  on  the  length  of  each  run.  Monitoring  was 
done  in  background  during  real-time  foreground  digitizing  at  a maximum  rate 


Spatial  Temperature 


Page 


8 


of  50  samples  per  second  per  channel.  Monitoring  consisted  of  dumping 
functions  of  various  seniors  to  strip  charts,  numeric  display  or  a printer. 
The  functions  generally  consisted  of  scaled  temperature  values  or  of 
differences  between  sensors.  The  selection  of  display  functions  and  sensor 
selection  was  completely  independent  of  foreground  digitizing  and  is  best 
categorized  as  a background  algebraic  language. 

DATA  ANALYSIS 


Introduct ion 

In  an  earlier  paper  ( Tillman,  J.  F-- , 1972  ) the  o /T*  relation  as  a 
function  of  Z/L  was  examined  with  the  aid  of  data  from  several  field 
experiments.  The  formulation  was  extended  to  cover  both  forced  and  free 
convection  as  opposed  to  the  prior  formulations  (Wyngaard,  et  al  . 1971) 
which  incorporated  only  free  convection  during  unstable  conditions.  Another 
limitation  of  these  results  is  that  they  include  only  "dry  convection",  i.e. 
those  cases  in  which  water  vapor  has  no  significant  effect  on  either  the 
heat  flux  or  the  stability.  In  the  present  paper,  the  results  are  extended 
to  cover  those  cases  in  which  the  latent  heat  flux  is  a major  component  of 
the  total  surface  to  atmosphere  heat  flux  and  in  which  water  vapor  is  a 
major  factor  in  determining  the  buoyancy  of  the  unstable  planetary  boundary 
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layer.  The  final  limitation  of  the  prior  efforts  is  that  the  height  of  the 
inversion  was  not  measured.  Consequently,  the  effect  of  7./?.^  on  the 
formulation  could  not  be  tested. 

Associated  with  the  more  complete  results  for  the  0-j/T*  versus  Z/L  relation 
are  the  surface  buoyancy  flux  results.  With  the  new  data  sets,  the  effects 
of  water  vapor  and  of  Z/Z^  on  the  calculated  buoyancy  flux  have  been 
estimated.  From  tills  analysis,  it  is  clear  that  the  surface  to  atmosphere 
total  buoyancy  flux  can  be  estimated  by  the  simple  measurement  of  the  mean 
and  standard  deviation  of  the  appropriate  temperature  throught  a major 
portion  of  the  planetary  boundary  layer. 

The  uniformity  over  space  and  time  of  several  temperature  temperature 
statistics  were  analyzed  using  the  array  sensors  during  the  period  of  from 
1500  to  1700  on  9/28/79.  This  was  a period  with  steady  light  winds  from  a 
favorable  direction  and  during  which  the  quiescent  temperature  intervals 
were  common.  Although  no  conaprisons  with  equivalent  measurements  from 
sensors  located  above  the  surface  layer  can  be  made,  the  observations  over 
time  and  space  are  used  to  inffer  some  characteristics  of  these  parcels  and 
t heir  origins . 
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The  definitions  and  symbols  for  this  paper  are  as  follows 


Z « height  of  observation 


Z^  = height  of  the  first  inversion 


T = average  temperature  over  time  or  space. 


When  latent  heat  flux  <<  sensible  heat  flux, 


temperature  is  used  in  all  statistics  while 


in  all  other  cases,  (unless  otherwise  specified) 


the  virtual  temperature  is  used  in  place  of 


temperature. 


T'  » deviation  of  temperature  from  the  average 


o =■  standard  deviation  of  temperature  over  time  or  space 


T*  " -H0/  p cpu*  " -w  T / u* 
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L = Obukhov  length 

kgZH0 

Z/L  “ Obukliov  stability  parameter 

I'q  = Surface-atmospheric  buoyancy  flux 

Hj,  = Vertical  buoyancy  flux  at  height  Z 
- P cpw  T 

k = von  Karman's  constant 
= 0-35 

Cp  = specific  heat  at  constant  pressure 

P = air  density 


u* 


* 

-u  w 
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p/T*  in  the  Planetary  Boundary  Layer 

In  the  prior  work  of  Tillman  ( 1972  ),  several  factors  in  the  °.j./T*  relation 
could  not  be  Investigated  due  to  limitations  of  the  data.  Since  this 
relation  is  the  basis  of  the  indirect  computation  of  buoyancy  flux,  it  is 
desirable  that  it's  range  of  application  be  as  wide  as  possible  and  that  its 
limitations  be  well  defined.  The  factors  that  were  either  poorly  or  not  at 
all  defined  were- 


The  need  for  a stability  correction  for  small  values  of  Z/L  where  both 
buoyant  and  mechanically  driven  turbulence  contribute  to  the  fluxes  of 
heat,  momentum  and  moisture , 

The  effect  of  moisture  on  the  buoyant  processes. 

> The  effect  of  the  height  of  observation  as  a function  of  the  inversion 

height,  7,/Z^.  and 

The  effect  of  the  actual  heat  flux  on  the  measurement  errors. 


To  understand  the  importance  of  the  parameters  not  previously  considered,  it 
is  appropriate  to  start  with  the  relation  between  °] /T*  and  Z/L  as  shown  in 
Figure  2,  incorporating  both  the  old  and  new  data  sets.  ( Note  that  the 
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height  range  O-CKZ/Z^O.S  is  presented  while  the  complete  layer  below  the 
inversion  0-0<7./Zi<1.0  is  presented  in  figure  3.  ) For  the  AMTEX , ROMEX . 

GATF  .Great  Lakes  and  Puerto  Rico  experiments,  the  data  were  obtained  from 
flights  using  National  Center  for  Atmospheric  Research  Aircraft  ( Hiester, 
1977,  Pennell,  W.  T.  and  M.  A-  LeMone.  1974-  These  data  are  hereafter 
referred  to  as  "W"  to  indicate  that  they  were  obtained  over  the  water.)  For 
the  Minnesota  experiment,  data  were  obtained  from  tower  instrumentation  up 
to  32  meters  and  from  balloon  instrumentation  up  to  the  inversion  height 
(Izumi  and  Caughey , 1976).  The  Kansas,  Round  Hill  and  Utah  experiments, 

KRU.  provided  the  data  for  the  earlier  results  ( Tillman,  1972,  Monji,  1972) 

Several  important  factors  are  evident  upon  close  inspection  of  this  figure. 

-1  /3 

The  first  is  that  the  KRU  data  closely  follow  a (7./L)  relation  which 

is  characteristic  of  free  convection  similarity  scaling.  The  only 
significant  deviation  is  at  the  small  values  of  7./L  where  both  forced  and 
free  convection  control  turbulent  interchange.  It  is  for  this  stability 
region  that  the  function: 

1)  °t/T*— 0-95(0. 0549-Z/L)"1/3 

was  formulated  ( Tillman,  1972  )•  The  Minnesota  data  were  added  to  verify 
the  relation  over  the  stability  range  previously  encountered  and  to  increase 
the  data  in  the  highly  unstable  region.  Since  this  is  the  region  where  free 
convection  similarity  scaling  and  the  -1/3  power  law  dependence  of  T*  on 
Z/L  should  be  most  accurate,  it  is  somewhat  surprising  to  find  the  large 
deviations  at  the  high  values  of  7./L.  The  first  explanation  is  to  assume 
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that  the  balloon  measurments  are  responsible  for  the  deviations,  where  in 
many  cases.  7.1  Z^  >>0.1.  As  the  convective  elements  loose  their  buoyancy 
excess  compared  to  the  surrounding  environment  due  to  entrainment,  inertial 
forces  become  dominant  and  "free  convective"  scaling  might  not  be  expected 
to  apply.  Partially  for  this  reason,  the  W data,  covering  the  range 
0 . 0<Z/ 7.  j < 1 . 0 , were  added  to  this  set  for  analysis  and  comparison.  Note  that 
even  though  the  W data  plotted  in  Figure  2 cover  the  range  0 . (KZ/Z^O . fl  . 
they  scatter  rather  uniformly  around  the  -1/3  power  law  relation  for  large 
~7./V.  rather  than  falling  below  the  free  convection  line  as  do  the  Minnesota 
values.  To  further  explore  the  Z/Z^  dependence,  the  same  data  are  plotted 
in  Figure  3.  stratified  by  Z/Z^,  except  that  the  points  now  cover  the 
complete  planetary  boundary  layer,  i.e.,  0 -0<Z/Z1<1 .0.  Note  that  with  the 
exception  of  the  Minnesota  data,  seven  points  in  the  range  O.R<Z/Zj <1 .0,  and 
a few  surface  based  points,  the  majority  of  the  planetary  boundary  layer  is 
well  characterized  by  the  "free  convection"  relation.  Consequently, 

one  is  led  to  conclude  that  the  deviations  of  the  Minnesota  data  at  high 
values  of  Z/L  from  the  "free  convection”  relation  can  not  readily  be 
attributed  to  Z/Z^  dependence  and  that  they  are  probably  due  to  instrumental 
or  meteorological  factors  unique  to  this  experiment.  It  is  fortunate  that 
the  -1/3  power  law  relation  holds  to  such  large  values  of  7.17,  , as  this 
allows  surface  buoyancy  flux  determination  throught  a major  portion  of  the 
boundary  layer  with  great  simplicity.  This  will  be  covered  in  the  next 
section. 


i 
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Oata  from  aircraft  measurements  over  the  water.  (Figure  4).  illustrate  the 
effect  of  water  vapor  on  the  °T/T<i  relation-  The  scatter  of  the  'dry' 
calculation  indicates  the  importance  of  water  vapor  in  the  majority  of  these 
data.  (The  error  is  so  large  for  the  dry  assumption  that  some  points  exceed 
the  upper  plot  boundary.)  That  the  other  experients  produce  a °t/T* 
relation  with  low  scatter  when  water  vapor  has  been  ignored,  must  be  due  to 
the  small  differences  between  statistics  of  temperature  and  virtual 
temperature.  Examination  of  data  from  the  Australian  experiment  (Shaw  1978) 
at  a site  with  a stubble  grass  cover  indicates  that  the  difference  between 
°T  an^  °TV  at  most  a few  percent-  However,  at  many  land  sites  and 
seasons,  the  latent/sensible  heat  flux  ratio  is  much  greater  than  one  (Munn, 
1966)  and  would  thereby  produce  significant  errors  in  any  °t/T*  versus  Z/L 
plots  derived  without  including  the  effects  of  water  vapor. 


Buoyancy  Flux 

A major  practical  reason  for  exploring  the  °-j /T*  relation  in  great  detail  is 
to  evaluate  the  possible  simplifications  and  the  limitations  when  using  it 
to  determine  the  buoyancy  flux  in  the  general  case  and  the  sensible  heat 
flux  in  the  case  of  dry  convection.  For  the  tower  based  surface  layer  data 
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previously  analyzed,  a stability  correction  for  the  stability  range  Z/L<0.2 
provided  increased  accuracy  in  the  heat  flux  calculations • (Tillman,  1972) 
However,  this  is  undesirable  in  that  the  heat  flux  can  no  longer  be 
determined  exclusively  from  constants,  pressure,  temperature  and  a, alone 
and  either  a direct  or  an  indirect  measurement  of  stability  is  required. 
Direct  measurement  of  stability  is  difficult  in  that  both  heat  flux  and 
shear  stress  measurements  are  required.  The  shear  stress  measurements  are 
noisy,  even  near  the  surface,  in  that  for  averaging  periods  of  less  than  an 
hour,  the  values  often  fluctuate  by  several  tens  of  percent.  As  one  departs 
from  the  proximity  of  the  surface,  the  scale  sizes  responsible  for  a major 
portion  of  the  stress  grow  larger  and  longer  averaging  times  are  required  to 
obtain  estimates  with  similar  accuracy.  In  addition  to  this  problem,  there 
are  a number  of  instrumental  problems  associated  with  the  stress 
measurement.  Indirect  determinations  can  be  made  using  relations  such  as 
those  relating  the  profiles  of  wind  and  temperature  to  the  fluxes  (Rusinger, 
1973)  or  such  as  the  skewness-Z/L  relation  (Tillman,  1972).  The  former 
require  several  levels  of  measurements,  although  they  seem  to  produce  stable 
estimate;  over  smooth  terrain-  The  latter  may  be  as  variable  as  the  direct 
measurements,  but  the  data  are  insufficient  to  verify  this  point. 
Consequently,  it  is  valuable  to  determine  the  accuracy  that  can  be  obtained 
by  ignoring  the  stability  correction. 

examining  the  complete  data  set  of  Figure  2,  it  is  evident  that  the  free 
convective  -1/3  power  law  fits  fairly  well  for  values  of  -Z/L>0.1.  To 
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determine  the  effect  of  ignoring  the  stability  correction,  first  consider 
the  comparison  using  this  correction,  important  for  the  small  values  of  Z/L. 
Figure  5 presents  the  ratio  of  the  calculated/measured  buoyancy  flux,  where 
the  calculated  complete  unstable  buoyancy  flux  relation  is  given  by: 


*,  l / 2 

2)  H/p  cp  = [(o  /Cl)J(kgZ/T)(C2-Z/L)/(-Z/L)] 


(For  this  calculation,  the  Obukhov  length,  L,  is  obtained  from  the 
measurement  of  it's  constituent  variables  while  in  the  earlier  work  by 
Tillman  (1972),  it  was  obtained  solely  from  temperature  statistics,  i.e. 
the  skewness  of  temperature-  Z/L  relation,)  Note  that  the  buoyancy  flux  is 
determined  through  the  complete  stability  range  and  most  of  the  planetary 
boundary  layer  without  major  bias.  If  the  Minnesota  data  are  excluded,  the 
results  are  even  better  and  the  majority  of  the  calculations  fall  within  -25 
to  +30  % of  the  measured  value.  Figure  6 presents  the  same  data  stratified 
by  Z/Zj.  The  results  indicate  an  underestimation  of  the  buoyancy  flux  for 
high  values  of  Z/L  if  the  Minnesota  data  are  included  whereas  the  relation 


does  not  appear  biased  if  these  data  are  excluded.  There  may  be  an  actual 


deviation  from  the  free  convection  relation  for  the  higher  values  of  7./L  but 
these  are  generally  also  the  cases  where  Z/Zi>>0-1  and  the  inertial  forces 
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the  complete  unstable  relation  is  shown  in  Figure  7.  There  appears  to  be  a 
tendency  for  the  error  to  decrease  as  the  magnitude  of  the  flux  increases; 
this  is  especially  true  if  the  Minnesota  data  are  eliminated  from  the 
analysis.  This  is  not  surprising  in  that  the  errors  in  measuring  the 
temperature,  and  the  buoyancy  flux,  should  increase  as  the  flux  decreases. 
The  plot  of  the  0,^,/T*  relation  of  Figure  4 also  illustrates  the  fact  that 
ignoring  the  effect  of  water  vapor  can  produce  major  errors  in  the  heat  flux 
determinations.  The  errors  assuming  only  "dry"  convection  are  large  enough 
for  the  over  water  experiments  to  completely  invalidate  the  relation. 

It  should  be  remembered  that  even  though  the  temperature  statistics  are 
obtained  at  a height  where  the  buoyancy  flux  is  a small  fraction  of  the 
surface  flux,  it  is  the  surface  flux  that  is  calculated  (7,*=0)  from  the 
temperature  statistics  at  height  Z.  If  the  buoyancy  flux  at  the 
observational  height  is  desired,  it  must  be  obtained  either  by  direct 
measurement  or  by  extrapolating  the  surface  value  to  the  height  of 
observation.  The  accuracy  of  this  extrapolation  depends  on  the  assumptions 
of  the  extrapolation  model.  Figure  8 compares  the  calculated  versus 
measured  "local  complete  unstable  buoyancy  flux"  for  the  W data.  In  this 
case,  the  extrapolation  is  given  by  the  simple  linear  extrapolation; 


3)  Hz  - H0(  -1.25  Z/Zi  + 1 ) 

where  the  buoyant  flux  goes  to  zero  at  Z/Zj  - 0-P.  (The  linear  decrease  of 
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buoyancy  flux  with  increasing  Z/Zj  is  implicit  in  the  assumptions  of  free 
convective  similarity  theory.)  As  expected,  the  calculated  local  values  of 
buoyancy  flux  are  less  accurate  than  the  surface  values,  the  error  becomming 
worse  with  increasing  Z/Zi-  An  important  fact  to  remember  is  that  in  the 
study  of  most  boundary  layer  processes,  eg.  by  modeling,  the  surface 
buoyancy  flux,  rather  than  the  local  buoyancy  flux,  is  a required  boundary 
parameter.  The  local  value  by  itself  is  not  particularly  useful  without 
either  a measurement  at  another  height  or  a knowlege  of  Z . in  addition  to 
the  height  of  measurement. 


The  "free  convection"  approximation  to  the  buoyancy  flux  provides  a good 
estimate  during  unstable  conditions,  when  the  flux  is  largest,  by  the  simple 
measurement  of  a temperature  statistic.  This  is  true  even  though  the 
temperature  measurement  is  made  well  above  the  surface  layer,  in  which  the 
fluxes  are  constant  with  height.  The  effect  of  ignoring  the  stability 
correction  to  the  buoyancy  flux  estimation  can  be  evaluated  by  considering 
the  analysis  of  Figure  8,  where  the  calculated  buoyancy  flux  is  given  by: 


4)  H/RHOcp  = [(uT/C1)3kgz/T) 


The  difference  between  the  complete  unstable  and  the  free  convection 
relations  is  less  than  20%  for  -Z/L>0.1.  It  departs  significantly  as 
neutral  conditions  are  approached,  i.e.  as  the  surface  is  approached  for 
fixed  buoyancy  flux  and  shear  stress  or  as  the  buoyancy  flux  decreases  or 
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the  shear  stress  increases  at  a fixed  height  during  differing  conditions. 

If  the  stability  range  is  restricted  to  -Z/L>0-1  and  the  questionable 
Minnesota  data  points  are  excluded,  a majority  of  the  estimates  are  within 
25%of  exact  agreement.  Examining  the  points  in  Figure  5 that  correspond  to 
the  values  of  Z/Z^>0.2  and  the  stability  range  0 . 1<-Z/I.< 10. 0 , the  majority 
of  the  largest  errors  are  from  the  Great  Lakes  and  Puerto  Rico  experiments- 
the  tower  data  are  well  behaved.  In  the  case  of  the  Great  Lakes  data,  the 
scatter  probably  is  due  to  inhomogeneity  and  non-stationarity  as  the 
experiment  was  conducted  to  study  the  effect  of  air  mass  modification  over 
the  Great  Lakes  during  frontal  conditions.  The  Puerto  Rico  scatter  is 
mainly  due  to  a faulty  air  temperature  sensor  and  the  fact  that  the  data 
agree  to  any  reasonable  extent  with  the  other  experiments  is  due  to  the 
dominance  of  the  water  vapor  contribution  to  Tv  and  the  buoyancy  flux. 
(Pennell,  private  communication.  1979) 


An  important  application  of  this  type  of  technique  is  the  measurement  of 
buoyancy  fluxes  where  the  instrumentation  packages  do  not  allow  the  cost  or 
complexity  of  direct  flux  measurements  or  their  inference  through  the 
flux-gradient  relations.  An  example  of  such  an  application,  is  the  Viking 
Mars  Lander  Meteorology  Instrument  where  wind  speed,  wind  direction, 
temperature  and  pressure  were  measured  at  one  height.  (Chamberlain,  et  al.. 
1976)  A comparison  of  fluxes  estimated  from  temperature  statistics,  using  a 
smaller  sample  set  than  optimum  and  using  the  free  convection  approximation, 
to  another  indirect  method  (Sutton,  Leovy  and  Tillman,  1978).  indicates 
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agreement  to  within  50%.  (As  the  latent  heat  flux  is  much  smaller  than  the 
sensible  heat  flux,  except  where  CO2  sublimation  is  taking  place,  the 
buoyancy  and  sensible  heat  fluxes  are  identical  in  the  absence  of 
sublimation.)  If  the  stability  dependance  and  other  factors  are  included, 
the  agreement  improves.  The  lack  of  better  agreement  between  the  two 
indirect  techniques  probably  is  due  to  the  limitations  of  the  current 
analysis  of  temperature  statistics  and  the  dissimilarities  of  the  data  sets 
analyzed  rather  than  inherent  disagreement  of  the  techniques. 

The  accuracy  of  the  complete  calculated/measured  buoyancy  flux  relation, 

(and  implicitly  the  o.j,/T*  relation)  is  impressive  if  it  is  realized  that  the 
data  were  taken  over  a 17  year  period,  by  numerous  investigators,  with 
differing  instrumentation,  observational  platforms  and  analysis  techniques, 
over  a height  range  of  0 -0<Z/Zj<0 . 8 , under  stationary  and  non-stationary 
conditions  as  well  as  over  homogeneous  and  inhomogeneous  terrain.  Another 
fact  to  remember  is  that  the  errors  in  measuring  stress  and  consequently  Z/l, 
and  T*  may  be  a major  factor  contributing  to  the  scatter.  Additional 
experiments  with  more  precise  measurements  of  heat  and  momentum  fluxes  and 
more  attention  to  the  other  factors  will  be  required  to  separate  the  effects 
of  the  various  parameters  and  improve  the  confidence  in  the  basic  relation. 
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Spatial  Characteristics  of  Convective  Elements 

One  of  the  interesting  characteristics  of  fast  temperature  measurements  in 
the  unstable  boundary  layer  is  the  periods  of  almost  vanishing  temperature 
fluctuations  when  the  surface  layer  is  super  adiabatic  or  even  dramatically 
auto-convective.  Even  though  the  lapse  rate  can  reach  almost  100  degrees 
Celsius  per  meter  in  the  lowest  meter,  (Ryan  and  Carroll,  1970)  periods  of  a 
half  to  one  minute  are  encountered  during  when  the  fluctuations  are  on  the 
order  of  a few  hundredths  of  a degree  or  less.  Parcels  with  sizes  such  as 
indicated  by  the  time  scale  must  originate  well  above  the  surface  layer  and 
be  advected  downward  into  the  surface  layer  by  the  large  scale  convection. 

It  is  suggested  that  the  large  scale  convection  is  responsible  for  a 
significant  portion  of  the  organization  of  the  smaller  elements  of  the 
convective  boundary  layer.  Although  the  demonstration  of  this  hypothesis  is 
beyond  this  analysis,  some  Insight  into  the  effect  of  large  scale  elements 
, on  surface  layer  can  be  obtained  by  considering  the  uniformity  over  time 

and  space  of  temperature  during  the  quiescent  intervals. 

During  stationary  conditions,  the  absoulte  temperature  within  the  quiescent 
intervals  repeatedly  reaches  the  same  minimum,  to  within  a few  hundredths  of 
a degree,  when  the  temperature  Is  devoid  of  fluctuations  greater  than  the 
order  of  0.01  degrees  Celsius.  If  a parcel  were  advected  from  a height  of 
several  tenths  of  the  inversion  height  and  the  height  is  on  the  order  of  a 
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kilometer  or  greater,  it  might  be  expected  to  have  a uniform  potential 
temperature  for  one  or  more  of  the  following  reasons: 

At  it's  source  height,  the  vertical  and  horizontal  gradients  would  be  small 
and  slowly  varying  in  time  so  that  parcels  advected  to  the  surface  sould 
initially  have  small  fluctuations. 

Any  fluctuation  would  be  decreased  during  the  downward  advection  due  to 
mechanical  mixing  except  at  the  boundaries  where  entrainment  is  occurring. 

If  the  statistics  of  Table  1 are  considered,  the  minimum  temperatures  of  the 
4 meter  sensors  is  quite  uniform  over  the  array  during  the  15  minute  period, 
varying  by  no  more  than  0.05  degrees.  However,  the  maximum  temperature 
varies  by  more  than  0.25  degrees  as  the  elements  are  quite  small  since  they 
originate  near  the  surface.  Table  2,  at  1719  shortly  after  the  transition 
to  stable  at  the  surface,  indicates  a reversal  of  the  roles  with  regard  to 
quiescent  temperature-  Since  the  surface  has  begun  to  cool  radiatively,  the 
potential  temperature  increases  with  height  and  the  parcels  advected  into 
the  surface  layer  are  now  warmer  than  their  surroundings  during  the 
quiescent  interval.  During  this  period,  the  range  of  maximum  temperatures 
over  the  horizontal  array  is  0-108  degrees  while  the  range  of  the  minimum  is 
0-268  ut,,rees. 

These  results  are  better  illustrated  in  Figure  9 which  summarizes  these 
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characteristics  for  each  period.  The  important  characteristic  to  note  is 
that  the  variation  of  the  minimum  ( and  equivalently  quiescent  ) temperature 
is  small  for  the  unstable  intervals  while  the  variation  of  the  maximum  is 
several  times  as  large.  However,  when  the  stability  reverses,  during  the 
transitional  period  between  16-07  and  17-19.  the  roles  reverse  with  the 
maximum  temperature  being  more  consistlent  in  time  and  space  and  the  minimum 
less  so.  (The  Minnesota  experiment  is  not  unique  in  illustrating  the  time 
variation  of  the  quiescent  temperature  by  itself.  However,  it's  combination 
of  time  and  space  measurements  with  high  accuracy  and  fast  response  is 
unique  except  for  the  recent  experiment  at  the  Roulder  Atmospheric 
Observatory  - ) 


From  these  and  other  observations,  (l.’ebb,  1965)  it  is  suggested  that  if  the 
quiescent  temperature  in  the  surface  layer  has  a minimum  that  repeats  to  an 
accuracy  of  a few  hundredths  of  a degree  within  a 15  minute  interval,  and  if 
the  fluctuations  are  on  the  order  of  a few  hundredths  of  a degree  or  less, 
then  the  air  parcel  originated  well  above  the  surface  layer  and  it's 
conservative  properties  are  characterises  of  the  majority  of  the  well  mixed 
convective  boundary  layer.  To  the  extent  that  this  is  true,  then  the 
measurement  of  conservative  properties  such  as  potential  temperature, 
humidity  and  particulates  during  these  periods  can  be  used  to  characteize  a 
significant  portion  of  the  planetary  boundary  layer.  Ry  combining  these 
with  mean  measurements,  which  are  characteristic  of  the  surface  layer,  is 
may  be  possible  to  monitor  the  evlution  of  some  variables  both  in  and  above 
the  surface  layer  from  surface  layer  measurements- 
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Table  1 


Temperature  statistics  over  15  minutes  from  the  University  of  Washington  and 
the  AFCRL  sensors.  Temperatures  are  given  in  degrees  Celsius  for  the 
parameters.  "T"  indicates  a UW  sensor  located  on  the  AFCRL  tower  at  the 
same  height  as  an  AFCRL  quartz  thermometer  mean  temperature  sensor.  The 
height  is  given  in  meters.  "C"  indicates  a UW  sensor  at  the  center  of  the 
array  (Figure  1),  located  50  meters  to  the  west  of  the  AFCRL  tower,  and  the 
height  is  given  in  meters.  The  other  designations.  "W",  "E" . and  "N" 
indicate  the  distance  to  the  west,  east  and  north  respectively  for  the 
sensor  as  referred  to  the  center  of  the  array.  All  of  these  sensors  are 
located  at  a height  of  4 meters-  Sensors  are  grouped  in  three  groups,  with 
duplicate  entries  in  each  group  when  appropriate  for  comparison.  The  first 
group  consists  of  all  sensors  on  the  AFCRL  tower  and  the  mean  temperature  of 
the  UW  and  AFCRL  sensors  is  identical  by  definition  for  this  table.  For 
other  time  periods,  these  mean  temperatures  should  agree  closely  if  the 
calibrations  are  accurate.  The  next  group  of  sensors  are  the  four  at 
different  heights  on  the  mast  in  the  center  of  the  array.  The  final  group 
include  all  sensors  at  the  height  of  4 meters,  irrespective  of  their 
location.  The  sampling  interval  is  15  minutes  over  which  45,000 
measurements  were  made  for  each  UW  sensor.  The  MFAN , MI\.  MAX,  SICMA  and 
SKEW  are  the  statistics  for  the  given  15  minute  period  without  the  removal 
of  any  low  frequency  components  or  trends.  Sensors  that  are  unreliable  are 
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indicated  by  comments  and  a thin  horizontal  line  through  the  data.  The 
skewness  is  the  normalized  third  moment  around  the  mean.  (Tillman,  1972) 
This  period  is  at  the  end  of  a moderately  stationary  unstable  daytime 
convective  time  interval. 

Table  2 

The  same  functions  as  Table  1 except  at  1719  the  surface  layer  has  now 
become  stable- 
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FIGURE  CAPTIONS 


Figure  (1) 


Schematic  of  array  layout  for  Minnesota  field  program.  Single  sensors  at  a 
height  of  4 meters  are  indicated  by  a dot  while  locations  with  multiple 
sensors  are  indicated  by  an  X-  At  the  array  center,  sensors  are  located  at 
heights  of  0.5,  1,  2 and  4 meters  while  on  the  AFCRI,  tower,  UW  sensors  are 
located  at  heights  of  2,4.8,16  and  32  meters,  along  with  AFCRL  sensors.  The 
horizontal  spacing  is  0.64,  2-7,  7.1,  20  and  50  meters. 


Figure  (2) 


o,j./T*  (where  T is  replaced  by  Tv  whenever  water  vapor  contributes  to  the 
total  buoyancy)  function  of  Z/L  stratified  by  experiment  for  the  normalized 
height  range  0 .0<Z/Zi<0 . 8 . For  the  AMTF.X , ROMF.X,  CATE,  Great  Lakes  and 
Puerto  Rico  experiments,  where  water  vapor  is  a major  factor  in  determining 
the  buoyancy  flux,  the  virtual  temperature  is  used  in  place  of  temperature. 
The  vertical  buoyancy  flux,  w'Tv'.  is  used  in  the  computation  of  L and 
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Figure  (3) 


°T /T * function  of  Z/I,  stratified  by  Z/Zt  for  the  normalized  height  range 
0.0<Z/7.i<l  .0.  The  characters  and  their  associated  ranges  are: 

0.0<Z/Z1<0.2=1 . 0 • 2<Z/7-i<0 . 4=2  , 0. 4<Z/7-i<0 . 6=3  . 0 . 6<Z/Z1<0 . 8=4  , 

0-8<Z/Zi<l -0=5.  The  data  are  characterized  by  their  range  of  Z/Zt  and  the 
range  0-8<Z/Z1<l-0  excluded  from  2 , is  included  to  illustrate  the 
increasingly  large  errors  encountered  as  the  inversion  is  approached-  these 
points  are  excluded  from  all  subsequent  plots. 

Figure  (4) 


°t/T*  as  a function  of  Z/L  for  the  AMTFX , BOMKX , CATF..  Great  Lakes  and 
Puerto  P.ico  experiments  for  the  normalized  height  range  0 .0<7./Zi<0 . 8.  The 
points  indicated  by  'V'  utilize  the  virtual  temperature  while  those 
indicated  by  *D*  use  only  dry  temperature. 
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Figure  (5) 


The  ratio  of  the  calculated  to  measured  buoyancy  flux  at  the  surface  as  a 
function  of  -Z/L  by  experiment,  computed  using  the  complete  unstable 
relation  , 2),  for  the  calculated  buoyancy  flux.  The  normalized  range  of 
heights  are  0.0<Z/Z^<0.8  for  land  and  oceanic  data  from  tower  and  aircraft 
platforms  respectively.  The  Minnesota  data  includes  both  tower  and  balloon 
measurements-  Where  the  bouyancy  is  significantly  affected  by  water  vapor, 
virtual  potential  temperature  replaces  temperature  both  in  the  direct 
statistics  and  in  the  computation  of  the  parameters.  Note  the  apparent 
trend  of  the  Minnesota  data  towards  an  underestimation  of  the  buoyancy  flux 
for  large  values  of  Z/Z^.  due  to  the  balloon  measurements. 


Figure  (6) 


The  ratio  of  calculated  to  measured  buoyancy  flux  as  a function  of  Z/L 
stratified  by  Z/Z^.  The  points  are  the  same  as  those  of  Figure  5 . The 
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range  of  normalized  inversion  height  represented  by  each  character  is 
specified  in  the  insert.  As  no  values  of  Inversion  height  are  available  for 
the  Kansas,  Round  Hill  and  Utah  experiments  and  they  were  made  from  towers 
less  than  32  meters  high,  all  data  from  these  experiments  are  assumed  to 
fall  in  category  1. 

Figure  (7) 


The  ratio  of  calculated  to  measured  buoyancy  flux  as  a function  of  the 
actual  measured  buoyancy  flux  for  the  height  range  0.0<Z/Z,j<0 -8- 

Figure  (8) 


The  ratio  of  complete  unstable  calculated/measured  buoyancy  flux  as  a 
function  of  Z/Z^.  The  data  are  clasified  as  surface  values,  dots,  and  local 
values,  crosses.  Local  of  buoyancy  flux  are  obtained  by  extrapolating  the 
surface  values  values  (determined  from  local  values  of  temperature 
statistics)  to  the  height  of  measurement  using  the  linear  extrapolation. 
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given  by  3).  Tower  measurements  are  excluded  from  this  plot  as  Z/Z ± is 
generally  unknown  and  the  local  and  surface  values  are  essentially 
ident ical • 

Figure  (9) 

The  ratio  of  calculated/measured  buoyancy  flux  as  a function  of  Z/L 
stratified  as  a function  of  Z/Z^.  The  calculated  buoyancy  flux  uses  the 
free  convection"  approximation,  given  by  4)-  Exact  agreement  for  the 
complete  unstable  rela  ionship  is  indicated  by  the  thin  lower  curve.  Data 
include  the  height  range  0 -0<Z/Z^<0 • 8 . 

Figuref 10) 

Statistical  properties  of  temperature  over  time  and  the  array.  The  top 
curve  illustrates  the  Obukhov  stability  paranerer  for  each  15  minute  run. 
going  from  unstable  at  the  beginning  to  stable  during  the  last  period.  The 
stability  values  are  derived  from  the  profiles  of  wind  and  temperature,  via 
the  flux-  gradient  relations.  The  second  curve  gives  the  mean  wind  speed  at 
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4 meters  during  the  same  Inteval-  The  group  of  three  curves  In  the  lower 
figure  all  present  data  from  the  15  minute  statistics  over  the  horizontal 
array  at  the  4 meter  height.  The  upper  of  the  three  curves  illustrates  the 
behavoir  of  the  mean  temperature  over  the  array  during  each  15  minute 
period-  For  each  time,  there  is  a minimum  and  a maximum  temperature 
displayed  along  with  a standard  deviation.  The  standard  deviation  is 
indicated  by  the  heavy  bar  which  is  PM  one  standard  deviation  around  the 
mean.  In  this  case,  the  mean  is  the  average  of  all  mean  temperatures  and 
the  standard  deviation.  Min  and  Max  refer  to  statistics  of  the  mean 
temperature  of  all  sensors  during  the  given  15  minute  period.  The  reason 
for  the  equivalence  of  the  data  for  the  mean  temperature  at  15:07  is  that 
the  sensors  were  defined  to  have  the  same  mean  temperature  at  tills  time  for 
cross  calibration  with  AFCRL-  The  minimum  and  maximum  temperature  plots  are 
constructed  in  the  same  manner  as  that  of  the  mean,  ie . from  the  minimum 
and  maximum  temperatures  at  each  location  during  the  15  minute  interval. 
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